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Abstract: 

The energy dependence for the survival probability of large rapidity gaps (LRG) 
<| S | 2 >, is calculated in an Eikonal model assuming a Gaussian opacity. The parameters 
determining <| S | 2 > are evaluated directly from experimental data, without further 
recourse to models. We find that <| S | 2 > decreases with increasing energy, in line with 
recent results for LRG dijet production at the Tevatron. 
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General approach 



Large rapidity gaps (LRG) have been studied both at the Tevatron Hp], and at HERA ||. 
LRG are expected whenever we have a process where a colour singlet is exchanged in the 
t channel. A prolific amount of data on LRG has been accumulated for dijet production, 
which is the channel that we will concentrate on in this letter. 

Historically, both Bjorken M and Dokshitzer et al. p[, suggested utilizing rapidity 
gaps as a signature for Higgs production, in the W-W fusion process in hadron-hadron 
collisions. As pointed out ||, the fraction of events for which a LRG is expected, is the 
product of two factors (see Fig.l) 

f gap =<\S\ 2 >-F s (1) 

where <| S | 2 > denotes the survival probability, i.e. the fraction of events for which the 
spectator interactions do not fill the rapidity gap of interest. F s is the fraction of events 
due to t-channel singlet exchange. Bjorken |4j] by comparing the rate for the exchange of 
two gluons, in a colour singlet state, to the exchange of one gluon in a colour octet state, 
estimated F s to be about 15 %. Other authors M M ]7[] M have suggested alternative 
ways to estimate F s , with results that are not too different from that of [|J. 

In this note we only concern ourselves with evaluating the first term in eq.(l) i.e. 
the survival probability. Experimentalists use a more practical definition of f gap . At the 
Tevatron f gap for LRG is defined to be the ratio: 

cross section for dijet production with LRG 
gap inclusive cross section for dijet production 

The jets are required to have a transverse energy E%r > 12 GeV at DO, and E J T et > 
20 GeV at CDF. 

The survival probability is easily defined using the eikonal model in impact parameter 
space. Normalizing our amplitude so that 

do~ , „, . ,n 

i '<«■«> I" 

<T„=4ir/m/(s,0) (2) 
The amplitude in impact parameter space is then given by: 

a el (s,b) = ±- [ d 2 qe^f(s,t) (3) 

Z7T J 
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from which we have 



o tot = 2jd%Ima el {sM 

a el = J d 2 b\ a el (s,b) \ 2 (4) 

s-channel unitarity implies that | a(s, b) \ < 1, and when written in diagonalized form we 
have 

2Ima el (s, b) =| a el (s, b) \ 2 + G in (s, b) (5) 
where the inelastic cross section 

a in = / d 2 bG in (s,b) (6) 



s channel unitarity is most easily enforced in the eikonal approach where, assuming that 
a(s,b) is purely imaginary, we have 

a el = i[l - e-^' 2 } 

G in (s,b) = l-e- n ^ (7) 

where Q(s, b) is a real function. 

Hence P(s, b) = e~ Q ( s,b > is the probability that no inelastic interaction takes place at 
impact parameter b. 

It is expedient to assume that the opacity b) can be approximated by a Gaussian 

n{s,b) = u(s)T(s,b) 

where 

v{s) = tt(s,b = 0) (8) 

with the profile 
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The Gaussian input eq.(9) corresponds to an exponential behaviour in t-space |]I0| . In 
the Regge framework the "soft" radius is of the form 

R 2 s (s) = A R 2 + A a (O)ln(-) = 2 B el (s) (10) 

So 

where B e i(s) denotes the forward slope of the elastic cross section (i.e. % ~ e Bel( - sS)t ). 
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With the amplitude in this form, the integrations in eq.(4) can be done analytically, 
and we obtain the following closed expressions for 



a tot = 2 / d 2 b(l - e~^> 6 )/ 2 ) = 2rrR 2 (s)[ln(^) +C- Ei(-^-)} 



'IV. 



d 2 b{\ - e - n{s > b) ) = nR 2 (s)[ln(u(s)) + C- Ei(-v(s))) (12) 



°ei = cr tot - a in = nR 2 (s)[ln(^Y) + C + Ei{-v{s)) - 2Ei{-^-] (13) 



where 



—dt 

-oo t 

and C = 0.5773. 

Note that the ratio only depends on the factor v(s). 

Following H the survival probability (see eq.(l)) <| S | 2 > is defined as the normalized 
product of two quantities. The first is the convolution over the parton densities, of two 
interacting projectiles presenting the cross section for the hard parton-parton collisions 
(dijet production) under discussion. The second term is P(s,b) i.e. the probability that 
no inelastic interaction takes place in the rapidity interval of interest. 

In the eikonal formalism 

2 Jd 2 bT H (b)P(s,b) 
11 Jd 2 bT H (b) { ' 

where P(s,b) = e~ n( - s ' b \ and the "hard" profile 

™ = ^JWf^' (15) 

where Rh denotes the radius of interactions in the "hard" scattering process. 
We define 



■<•> - M (16) 



i.e. 

, . 7cR%(s) interaction area for soft collisions 
a ( s \ — 15 v ' = i i > i 

7rR 2 H (s) interaction area for hard collisions 
Evaluating eq.(14) yields 

where the incomplete gamma function 7(0, x) = Jq z a ~ 1 e~ z dz 



<\S\ 2 >= ^W^4^ (17) 
1 1 [u(s)] a ^ 
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Numerical evaluation of <| S | 2 > 

Previous estimations of the LRG survival probability @] || were derived from various 
models assumed to describe the "soft" and "hard" interactions. Our approach here is to 
determine the parameters on which <| S \ 2 > depends in the Eikonal Model, directly from 
the relevant experimental data without further recourse to models. 

The survival probability is sensitive to the value taken for the radius of the "hard" 
interaction Rh- We wish to evaluate this radius directly from the experimental data, 
namely, from the cross section of the double parton interaction measured by CDF at 



the Tevatron JTT|. Indeed, the CDF collaboration have measured the process of inclusive 
production of two "hard" jets with almost compensating transverse momenta in each 
pair, and with values of rapidity that are very similar. Such pairs can only be produced 
in double parton collisions, and their cross section can be calculated using the Mueller 
diagram given in Fig. 2. 



The double parton scattering cross section can be written in the form ( see Fig.2) [ |TT|j : 

a ind {2 jets) a ind (2jets) 
a DP = m , (18) 

2(X eff 

where the factor m is equal to 2 for different pairs of jets, and to 1 for identical pairs. 
The experimental value of a e ff = 14.5 ± 1.7 ± 2.3 mb. In Fig.2, a e ff is obtained by 
integrating over momentum transfer q, or over the impact parameter b, namely 

J- = / d 2 bT 2 H (b) = -Lp . (19) 

Comparing with the experimental value of o" e //, one obtains R 2 H = 5 — 7 GeV~ 2 . 

Recent data from HERA |12| indicate that the photo and electro production of J/ip 
is a "hard" (short distance) process. Data has <Jh{i + P — y J/ip + p) ~ W ' 8 while the 
analagous "soft" process a s {^ + p — * p + p) ~ W^ 22 . The forward slope -B e z(7* + P — ► 
J/ip + p) ~ 4 GeV~ 2 , while for the inelastic diffraction process 7* + p — > J ftp + X, 
Btnil* + P — ► J ftp + X) ~ 1.6 GeV~ 2 , both B e i(s) and B in (s) appear to be independent 
of Q 2 and W, as one would expect for a "hard" process. This experimental observation 
may suggest that the proton is a two radii object |I3| . For the present calculation we use 



an average radius which we estimate to be R 2 H = 8GeV~ 2 . This value is smaller than 
than has been assumed up to now || 0. This is a conservative choice, and it is possible 
that future experimental data will dictate an even smaller value of R 2 H = 5 — 7 GeV~ 2 , 
as the CDF data for the double parton cross section suggests. 
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The pp and pp channels have measurements of B e i(s) = (| i? 2 ) and over a large 



energy range [|TJJ]. The parameter a(s) = R ^ can thus be determined directly from 



experimental data. The second parameter which is required to evaluate <| S \ > is u{s) 
(see eq.(8)). Using experimental values from pp and pp scattering data [14| for the ratio 
we can find the corresponding value of u(s) at the appropriate value of s ( energy 
squared), using eqns.(ll) and (13). The dependence of v(s) as a function of the ratio 
for pp scattering is shown in fig. 3. 

Once we have determined z/(s) and a(s) (eq.(16)) , the survival probability can be 
determined from eq.(17). In Fig.4 we display <| S | 2 > as a function of z/(s) and a(s). In 
Fig. 5 we show a contour plot for <| S | 2 >, for various values of v(s) and a(s). 



Comparison with data 



Recently DO [[Uj have measured LRG in dijet production in pp scattering at center of 
mass energies of ^fs = 630 and 1800 GeV, and find f gap {\fs = 630GeU) = 1.6 ± 0.2% 
and fg ap (^s = 1800GeV) = 0.6 ± 0.1% Q i.e. 

R f g g P (VS = 630 GeV) 

" f gap (V~s = 1800 GeV) ~ ^ ± U ' b 

The ZEUS collaboration reported a value of f gap =7 ± 2 %, at a center of mass 
energy < W >~ 150 GeV in the photoproduction of dijets. 

The survival probability <| S | 2 > increases as we go to lower energies, which corre- 
sponds to lower values of z/(s) and a(s) (see fig.4). In fig.5 we have marked the values of 
u(s) and a(s) which correspond to \/s = 546 and 1800 GeV. We would like to stress that 
the input data used in calculating <| S \ 2 > was taken from p-p and p-p interactions. Our 
predictions for the LRG survival probability are: 

(GeV) <| S \ 2 > 

150 16.3 % 

630 12.8 ± 0.4 % 

1800 5.9 ± 0.4 % 



Hence the ratio of: 



< 


s 


> x /i=630 


< 


s 


>7i=1800 



2.2 ±0.2 



*Prelimenary results presented by J. Perkins |T(| at the LISHEP98 Workshop reported a new value 
of R = 3.4 ± 1.3 
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Conclusions 



Recent measurements at HERA 0] provide us, for the first time, with a value for the 
radius of the "hard" interactions i.e. about 8 GeV~ 2 . Assuming a one radius model for 
the proton, and a Gaussian form for the profile in impact parmater space, gives us an 
explicit expression for <| S \ 2 >. 

We note that the LRG survival probability <| S \ 2 > decreases with increasing energy 
(see Fig. 5). This is a consequence of the fact that experimental data has both 
and B e i(s) increasing with s. This trend is consistent with the observations of the DO 
collaboration Jl5|]Jl6|] who find 

MVZ = MOGeV) _ 34± 13 



f g a P {V^ = 1800GeV) 
Using exprimental data to directly fix the parameters v(s) and a(s) we find 

2.3 ±0.2 



< 




2 \ 

> v /J=546GeV 


< 


s 


> y/S=1800GeV 



which we can interpolate to values measured at Tevatron energies i.e. 

2.2 ±0.2 



< 


s 


I >Vi=630Gey 


< 


s 


2\ r- 



This is consistent with F s , the singlet to octet colour exchange contribution (see eq.(l)) 
being mildly energy dependent as suggested by Eboli et al. [0, or energy independent as 
proposed by Buchmuller and Hebecker ||. 

The main results of this paper, that the survival probability decreases with increasing 
energy, and the predicted rates, are in agreement with what is observed at the Tevatron. 
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f = 



o (LRG) 



9 a P G(INCL) 



= <S > 




2 




Figure 1: Pictorial definition of / 9ap , where P and G represent respectively, the exchange 
of a colour singlet and a colour octet 
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Figure 2: The Mueller diagram for the double parton interaction. . 
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